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A method for estimating the duration of isothermal sintering of a powdered body and its final porosity and 
shrinkage is considered. The method, based on a three-dimensional model in the form of a cubic packing of 
spherical particles, uses sintering physics equations, the physicochemical properties of materials, and particle 
sizes. The model allows for describing not only the initial state (coalescence of particles) and the final state 
(compression of individual pores), but also the intermediate consolidation stage according to the general mass 
transfer mechanism. The sintering parameters for a sample of container glass powder are calculated using the 
viscous flow mechanism. The absolute discrepancy between the estimated results and experimental data is 
equal to 2.0% in porosity and 1.5% in shrinkage. 


The quantitative description of the consolidation of pow¬ 
dered bodies in sintering is not only of theoretical interest but 
of great practical significance as well. The theoretical esti¬ 
mate of the duration of sintering powder preforms to a preset 
density makes it possible to determine the firing temperature 
and the interval of the sintered state for particular materials, 
as well as the effect of the phase composition and the struc¬ 
ture of the powder. Such estimate can be obtained using a 
calculation method based on the physical regularities of mass 
transfer, whose parameters are the physicochemical proper¬ 
ties of materials and the structural parameters of particles 
and their packing. 

The main principles of the physics of sintering proposed 
in the works of Ya. I. Frenkel’, G. S. Kuchinskii, B. Ya. Pi¬ 
nes, Ya. E. Geguzin, et al. are summarized in [1]. The kine¬ 
tics of sintering of various materials is described by uniform 
equations whose parameters are determined by the prevailing 
mechanism of mass transfer. For most mechanisms (diffu¬ 
sion across the volume of crystals, along the boundaries, 
along the surface, viscous flow, transfer via the gaseous 
phase, sintering under an external mechanical load) kinetic 
equations have been derived and experimentally validated 
and the conditions of their applicability have been identified. 
However, the use of these equations in practical calculations 
is limited, since they are based on unidimensional flat mo¬ 
dels and describe only the initial (coalescence of particles) 
and the final (compression of isolated pores) sintering stages. 
The intermediate stage of a decreasing volume of communi- 
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eating permeable pores in a contiguous solid-phase matrix 
has not yet been described. The known estimates [2] of iso¬ 
thermal sintering duration based on one-dimensional models 
differ from real values a thousand times. 

In the context of the phenomenological approach, sin¬ 
tering can be described as the process of consolidation of a 
real powder body within the entire porosity variation interval 
[3]. The extension of this method has led to the development 
of methods for the mathematical modeling of sintering [4]. 
The parameters in the phenomenological equations are em¬ 
pirical ones determined experimentally. Lately attempts have 
been made to correlate the physical and the phenomenologi¬ 
cal approaches by estimating the parameters of the pheno¬ 
menological equation based on the physicochemical proper¬ 
ties of the sintering material [5]. 

In applied research the theoretical conclusions are 
mainly used qualitatively. 

For estimating the total duration of the consolidation of a 
powder body using the equations of the physics of sintering, 
it is necessary to develop a unified physical model of the ma¬ 
terial structure at all stages of the process, including the in¬ 
termediate stages, taking into account the volumetric packing 
of particles. 

We have earlier [6] proposed a 3D model of a powder 
body, the sequence of the structure modification stages in 
sintering, and methods for calculating the duration of these 
stages based on kinetic equations of the physics of sintering. 
These methods were used to describe the diffusion sintering 
of corundum ceramics. 
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Fig. 1. Cubic packing of particles. 



Fig. 2. Face of the cell of a sample at the stage of sintering [the ver¬ 
tices of the square correspond to the centers of particles; dashed 
lines) interparticle boundaries]. 


The purpose of the present study is to evaluate the suit¬ 
ability of this method for describing viscous-flow sintering. 

Three-dimensional packing model and sequence of sin¬ 
tering stages. The following assumptions have been made: 

- all particles are spherical of radius R ; 

- the particles form an elementary cell that can be de¬ 
scribed as a regular polyhedron; the size of the cell is signifi¬ 
cantly smaller than the size of the body; 

- the coordination number (CN) of particles in sintering 
does not change; 

- long-range arrangement of particles is absent; i.e., the 
symmetric axes of the cell are oriented with respect to the 
symmetry axes of the body at different angles with equal 
probability; shrinkage is isotropic; 

- the body structure is homogeneous; 

- sintering in the entire porosity interval occurs by the 
same mechanism in isothermal conditions; 

- sintering proceeds by the most favorable way; i.e., 
such undesirable phenomena as localization of sintering, co¬ 
alescence of pores, etc. are absent. 

The simplest 3D model of the packing of spherical parti¬ 
cles is the cubic model (CN = 6) whose porosity (47.64%) 
agrees with the porosity of actual preforms made of narrow- 
fractionated powders (45 - 55%). 

Figure 1 shows the section of a cubic cell across its face 
in the initial state. Parting from the geometry of the model 


cell, the shrinkage along its edge (convergence of particle 
centers) is equal to 

A/// 0 = 1 - cos (arcsin ( x/R )), (1) 

where x is the radius of the contact bridge between the parti¬ 
cles (Fig. 2). 

The kinetics of the bridge growth is described by the ex¬ 
ponential equation of the following form: 

x n - Xq = (x - X o)f(A, a, R m , 7), (2) 

where A is the mass transfer kinetic parameter; a is the spe¬ 
cific surface energy; x is the time; T is the thermodynamic 
temperature; n and m are constants. 

The beginning of sintering (x 0 = 0) corresponds to point 
contact (x 0 = 0). 

Equation (2) in the case of a viscous flow has the follow¬ 
ing form [1]: 
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where rj is the viscosity of the material at the specified tem¬ 
perature. 

The estimated dilatometric (macroscopic) shrinkage of 
the sample AT/Z 0 and the convergence of the particles cen¬ 
ters in the cell A/// 0 are proportional: 
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In view of the assumption on the absence of the 
long-range order in the body structure, the coefficient K in 
the case of the simplest cubic packing is equal to the mean of 
the lengths of the cube edge and the cube diagonal, i.e., 

1(1 + a/3) » 1.366. 
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Within a single cell volume, the volume of the solid 
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phase is equal to — nR 3 . The porosity of the sample is related 
to shrinkage: 
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At the first state of sintering, particles coalesce, i.e., their 
centers converge (shrinkage) and the surface area of the 
bridges between the particles grows. 

According to expression (4), dilatometric shrinkage of 
5.00% in the cubic packing corresponds to shrinkage within 
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a two-particle system of 3.66%, x 1 /R = 0.26808 according to 
Eq. (1) and the porosity according to expression (5) is equal 
to Por { = 0.414 or 41.4%. After a transformation according 
to Eq. (3) it is possible to determine the duration of the first 
state t 1 . A shrinkage equal to 5% is accepted, since accord¬ 
ing to the calculation, the total sintering duration in the linear 
shrinkage interval A L/L 0 from 2 to 15% varies within the 
limits of ± 7% of the process duration in relation to the 
shrinkage value at (A L/L 0 ) x = 5%. 

By the moment t 1 the maximum surface curvature is reg¬ 
istered in the bridges between the particles and the sites 
around the interpore bridges (“the pore neck”). The curvature 
radius of the “pore neck” r' on the face of the cell with the 
edge length 2R [1 - (A/// 0 ) 1 ] connecting the central hollow 
spaces inside the neighboring cells (Fig. 2) is equal to 

r , = R V-(M/lo)l\ _ R (6) 

cos(ti / 4) 

At the second stage, the convergence of particle centers 
and the isolation of pores proceed independently. One needs 
to take into account the isolation of pores at this stage and 
their subsequent spheroidizing for transition from a system 
of contacting spheres of the solid phase to an ensemble of 
isolated pores in the solid matrix. The solid phase surface on 
the cell face has the maximum positive curvature, since the 
growing bridges of the neighboring particles are maximally 
close. The isolation of pores follows the mechanism of sur¬ 
face mass transfer. 

Since other mass transfer mechanisms, apart from vis¬ 
cous flow, are absent from glass powder sintering, this pro¬ 
cess is also described by Eq. (3). 

At the initial moment, a sphere of the radius r' can be in¬ 
scribed into the pore neck (Fig. 2). At the end of the process 
the entire face of the cell exists in the solid phase of closely 
packed particles. The interface surface in the adjacent sites is 
formed by conjugated convex and concave surfaces of the 
particles and the emerging bridge. In the center of the face, 
slightly lower and slightly higher than the plane shown in 
Fig. 3, the solid phase has a spherical surface, concave to¬ 
ward the face plane. 

Let us imagine that at the moment when the sealing of 
the pore neck ends, two spheres of the radius r' (Fig. 3) con¬ 
tact at a point in the center of the face. The spheres are hol¬ 
low. The process of sealing of the pore neck can be repre¬ 
sented as the coalescence of two hollow spheres in reverse 
time, i.e., at the initial conditions x' Q) =r' with t 0 = x' and the 

final conditions x' = 0 with x = 0. We replace the radius R in 
Eq. (3) by r'. The solution of this equation yields the value x'. 

By the end of sealing of the central pore neck x 2 = x', the 
radius of the bridge according to Eq. (3) is equal to x 2 , the 
estimated dilatometric shrinkage according to formula (4) is 
(AL/Z 0 ) 2 , and the porosity is Por 2 . 



Fig. 3. Halves of hypothetical hollow spheres in the center of the 
bottom face of the cubic cell [dot-and-dash lines) cell edges]. 


We assume all pores to be sealed after the sealing of the 
pore necks (at the moment x = x 1 + x 2 ). Actually at the speci¬ 
fied level of shrinkage and porosity some pores remain open. 

Let us estimate the part of the interpore bridges. Each 
pore contacts six neighboring ones. The ratio of the volumes 
of the adjacent segments F segm to the total volume of the 
pores is 


6V. 


segm 


Por 2 [2R(l -(Al/ / 0 ) 2 )] 


0.194. 


According to the calculation Por 2 « 30% and the volume 
part of the interpore bridges is 0.3 x 0.194 = 0.0582, i.e., is 
relatively low; therefore, we ignore the error in describing 
the structure of these sites. 

The structure of a sintering glass powder preform is cha¬ 
racterized by the absence of boundaries between particles, 
i.e., we assume that after the isolation of pores the individu¬ 
ality of particles is lost and the pores are dispersed in the 
continuous glass matrix. Therefore, when estimating the 
shrinkage and consolidation of the body at the stages of 
spheroidizing and pore compression, one should take into ac¬ 
count the increased system viscosity due to the inclusion of 
the second phase, namely, the pores. For this purpose we 
used the dependence of viscosity on the volumetric share of 
the dispersed phase C v [7]: 
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where r| 0 is the viscosity of the dispersion medium (material) 
at the specified temperatures; K c is a constant depending on 
the shape and material of the dispersed phase inclusions; and 
C v cr is the maximum content of the disperse phase in the sys¬ 
tem being deformed. 

Let us assume K c = 2.5 for isometric inclusions, 
C v = Por 2 , and C v cr = V2 ^ « 0.7405, which corresponds to 

the degree of filling of the densest face-centered cubic pack¬ 
ing of spheres. 
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Fig. 4. Section of a cell by plane (101) at the stage of spheroidizing 
of the pore (a) and at the state of compression of the isolated 
pore ( b ). 


The presence of the second phase inclusions does not af¬ 
fect the rate of changing of the interphase surface geometry 
but modifies the kinetics of shrinkage and consolidation at 
the later stages of sintering. Therefore, the duration of the 
pore insulation stage (x 2 = x') and the spheroidizing stage x 3 
are calculated using the value r| 0 in Eq. (3), which corre¬ 
sponds to the viscosity of pore-free glass: 
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At the third stage of sintering, isolated pores become 
spheroizided due to the different curvature of the convex sur¬ 
faces of particles directed inside the central pore and the con¬ 
cave surfaces in the area of the cell face center near the pore 
neck. 

At the third stage, porosity and shrinkage are calculated 
based on the same equations as at the preceding consolida¬ 
tion stages. 

The geometrical considerations (Fig. 4 a) yield the fol¬ 
lowing relations between the thickness of the interpore 
spaces x ", the distance between diagonally positioned parti¬ 
cles X, and the angle between the segment connecting the 
centers of the particle and the conjugacy point of the convex 
and concave surfaces of the interface cp: 

x " = R sin cp - [R (1 - sin cp) + X/2] -—; (9) 

coscp 

, — *(1— sin(p) + X/2 (1Q) 

coscp 


The maximum curvature is observed on the surface of 
the interpore space (cup) between the particles located on the 


lateral diagonal of the cell. Spheroidizing ends when its cur¬ 
vature radius r becomes equal to the pore radius r por : 
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The radius r is determined based on Eqs. (9) and (10) 
using formula (8) after the calculation of x 

The solution of Eq. (10) with respect to cp at r = r and 
X/2 = r' yields the value of the angle cp end (the moment of the 
end of spheroidizing). The substitution of cp end into expres¬ 
sion (6) makes it possible to find the value x " nd and the solu¬ 
tion of Eq. (8) for x £ = 0 and x " nd gives the duration of the 

pore spheroidizing stage x 3 . 

According to Eq. (3) with a correction for viscosity vari¬ 
ation based on formula (7), within this time the value x grows 
to x 3 , which, considering expressions (4) and (5) corre¬ 
sponds to the shrinkage (AT/Z 0 ) 3 and the porosity Por 3 . 
Then the radius of spherical pores r por sf is recalculated sub¬ 
stituting (A/// 0 ) 3 in formula (11). 

The structure of the sample at the beginning of the final 
stage of sintering represents a cubic packing of isolated 
spherical pores of the radius r por sf in a continuous matrix. At 
this stage the pores are compressed under the effect of capil¬ 
lary pressure P cap and external pressure, which are opposed 
by the gas pressure P gas in the pore. The compaction ends at 
the moment of reaching their equilibrium. The final shrink¬ 
age, porosity, and pore size for an insoluble gas depends on 
its quantity contained in isolated pores. 

The equilibrium state of spherical pores is characterized 
by the radius r end calculated by the Boyle - Marriott law, 
assuming the gas in the pore to be ideal (Fig. 4 b ). 

For a cell with a single isolated pore the following equa¬ 
lity holds: 


jz gas = 
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Vf* s =V 3 Por 3 , 

where V 3 = 8R 3 [1 - (A/// 0 ) 3 ] 3 and Por 3 are the volume of 
the cell and the porosity at the moment of the end of pore 
spheroidizing; P 3 is the gas pressure at the moment of isola¬ 
tion equal to the pressure in the furnace. 
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The effect of external pressure on the value r por end is ne¬ 
glected. In the state of equilibrium 
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The equilibrium pore radius is equal to 
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The duration t 4 of the final stage of compression of iso¬ 
lated pores from the radius r por to r por end is calculated based 
on the equation below, which is common for the cases of 
small and large P (compared with capillary pressure P cap ) [1]: 
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where P is the external pressure; r| 4 is the viscosity of the po¬ 
rous medium based on Eq. (7) at C v = Por 3 . 

The total shrinkage of the model cubic cell is calculated 
based on formula (11) by substituting the value r por end , and 
the final porosity is found based on Eq. (5). 

The total duration of sintering in isothermal conditions is 
equal to the sum of the durations of individual stages 
(t 1 + t 2 + t 3 + t 4 ). 

Calculation of sintering duration . The material of the 
samples is clear bottle glass of the following composition 
(wt.%): 70.56 Si0 2 , 0.32 A1 2 0 3 , < 0.06 Fe 2 0 3 , 9.39 CaO, 
3.32 MgO, 15.66 Na 2 0, and 0.64 S0 3 . 

We assume the samples to consist of spherical particles 
of the radius R = 10 pm, the glass is a Newtonian fluid [8], 
its sintering proceeds in isothermal conditions in air under a 
pressure of 10 5 Pa, and the air is insoluble in the material. 

The temperature dependence of glass viscosity [9] in the 
range of 10 L2 ° - 10 4 - 97 Pa • sec at the temperatures of 
800 - 1400°C is approximated by the following polynomial: 

t] 0 = -4.9999 x 10 -9 ^ 3 + 2.28691 x 10"V- 
3.7990 x 10-^ + 23.291. 

The specific surface energy (surface tension) of glass at a 
temperature of 1300°C calculated according to the additivity 
rule [9] is a 1300 = 0.318 J/m 2 , its temperature dependence 
a ft) was assumed to be linear: 

a (0 = 0.318-3.18 x 10“ 5 (^- 1300). 

At the temperature of 700°C this glass is characterized 
by rj 0 = 1.5 x 10 7 Pa • sec, a = 0.337 N/m 2 . 

The duration of the first stage of coalescence at the tem¬ 
perature of 700°C corresponding to the linear shrinkage of 
samples (A L/L 0 ) l = 5% and Por ] = 41.4% according to 


TABLE 1 


Estimated parameters of glass powder samples 


Temperature, 

°C 


with particle diameter, pm 


5 

10 

20 

30 


Duration of isothermal sintering, min 


400 

2.9 x 10 5 

4.3 x 10 5 

5.2 x 10 5 

5.5 x 10 5 

450 

2.6 x 10 4 

3.7 x 10 4 

4.5 x 10 4 

4.8 x 10 4 

500 

2.7 x 10 3 

3.9 x 10 3 

4.7 x 10 3 

4.9 x 10 3 

550 

3.3 x 10 2 

4.8 x 10 2 

5.8 x 10 2 

6.1 x 10 2 

600 

49.0 

71.0 

86.0 

90.0 

650 

8.5 

12.0 

15.0 

16.0 

700 

1.7 

2.5 

3.0 

3.1 

750 

0.40 

0.57 

0.69 

0.72 


Shrinkage A L/L Q 



400 

25.92 

24.86 

21.98 

18.46 

750 

25.86 

24.78 

21.77 

18.08 



Porosity, % 



400 

1.6 

4.3 

11.4 

19.1 

750 

1.7 

4.5 

11.9 

19.8 

Eqs. (1), (3), and (4) is 

equal to t 1 = 

0.36 min. The pore neck 


radius r' = 3.63 x 10 6 m. 

The duration of the second stage (isolation of pores) is 
equal to u 2 = x' = 0.66 min, the linear shrinkage of the sam¬ 
ple is 14.7%, and the porosity Por 2 = 26.4%. 

The duration of the stage of spheroidizing t 3 = 1.53 min. 
Taking into account formula (7), according to Eqs. (1), (3), 
(4), and (5) the end of the third stage corresponds to 
(AL/L q ) 3 = 20.5%; Por 3 = 14.8%. After the recalculation, 
the radius of the spherical pores r por sf = 7.10 x 10" 6 m. 

At the moment of the end of sintering the equilibrium ra¬ 
dius calculated based on formula (13) r por end is equal to 
5.11 x 10“ 6 m. Taking into account expression (7) from 
Eq. (14), we obtain the duration of the pore compression 
stage t 4 = 0.42 min. 

Thus, according to the proposed calculation method, the 
total duration of isothermal sintering of a glass powder sam¬ 
ple via the viscous flow mechanism at a temperature of 
700°C is equal to about 3.0 min, the final linear shrinkage of 
the sample is 21.8%, and the porosity 11.8%. 

Let us estimate the effect of the dispersion of powder on 
the estimated values of the final shrinkage and porosity, as 
well as the sintering duration at different temperatures 
(Table 1). 

The pore volume distributed between the particles in 
sintering is related to the size of this particle by relations (5), 
(11), and (13). The weight of the gas contained in a pore at 
the end of the third stage depends on the size of the cell, i.e., 
on the size of the initial particles. The amount of gas in iso¬ 
lated pores is constant; therefore, the smaller the initial parti¬ 
cle size, the smaller the volume of captured gas, the equilib¬ 
rium pore radius, and the porosity of sintered material. 

The viscosity of sintering material determines the dura¬ 
tion of sintering. The effect of temperature on final shrinkage 
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Fig. 5. Temperature dependence of the final shrinkage of samples (a) 
and shrinkage kinetics of glass samples at 700°C ( b)\ dots) experi¬ 
mental data; solid curves) calculated values. 


and porosity is determined by the weak temperature depen¬ 
dence of the specific surface energy of glass and is insignifi¬ 
cant compared to the effect of the particle size. An increase 
in the pore size from 5 to 30 pm raises the duration of sin¬ 
tering by less than two times. The change in the shrinkage 
and porosity of sample within the firing temperature interval 
of 400 - 750°C does not exceed the experimental error in de¬ 
termining these parameters. 

Experimental validation of the method. The experiment 
was performed on glass powder of the specified composition 
ground by wet milling in a ball mill until completely passing 
through a sieve with a cell size of 40 pm. The fraction of 
8-12 pm with a mean particle radius of 10.4 pm was sepa¬ 
rated by sedimentation. Disk-shaped samples of diameter 
16 mm and height 2-3 mm were molded at a pressure of 
20 MPa. 

The porosity of the preform was 44.3%, which is close to 
the density of the simplest cubic packing. 

The method of consecutive firings in an electric furnace 
in air (Fig. 5a) was used to investigate the temperature de¬ 
pendence of shrinkage. The heating rate was 2.5 K/min. The 
exposure duration at the maximum temperature was 6 min 
including the estimated duration of heating the sample in the 


mold, which at 500-800°C was equal to 2.8-3.2 min. 
Cooling was performed together with the furnace. 

The kinetics of sintering was studied at the temperature 
of 700°C (Fig. 5b). The process mainly ends within a time 
period close to the theoretical estimate, taking into account 
the preform heating period (x « 6 min). When the exposure 
was extended to 120 min, the shrinkage of the samples grew 
by no more than 0.6%. The difference between the estimated 
ad experimental shrinkage is not more than 1.5%. 

The difference between the experimental isothermal 
sintering duration and the duration calculated using the 
above method is no more than tens of percent, whereas cal¬ 
culations based on the one-dimensional model yield a dura¬ 
tion value that is a thousand times longer. 

The experimental values of porosity (10.3%) and shrink¬ 
age (23.1%) of samples after firing with 6 min exposure at 
700°C agree satisfactorily with the calculated values of 
12.3% and 21.6%, respectively. The deviations are related to 
the fact that the porosity of the actual preforms (44.3%) is 
somewhat lower than the porosity of the model cubic cell 
(47.6%). 

Thus, the proposed model and the calculation method 
based on this model adequately describe sintering of glass 
powder samples by the viscous flow mechanism in isother¬ 
mal conditions. The calculation makes it possible to estimate 
the duration of isothermal sintering, the final shrinkage, and 
the porosity with an accuracy sufficient for practical use. 
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